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Abstract 

Dynamical focusing of ensembles of neutral particles in energy and configuration space has been 
demonstrated recently [C. Petri et al. 2010, Phys. Rev. E (R) 82, 035204] using time-dependent 
elliptical billiards. The interplay of nonlinearity, dissipation and driving yields the occurrence of 
attractors in the phase space of the billiard. Here we show that dissipative oval billiards with slowly 
oscillating elliptical scatterers in the interior allow for a dynamical focusing on simple periodic 
trajectories with close to perfect efficiency. This setup should be more amenable to corresponding 
experiments which are briefly discussed. 
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Introduction Exhibiting typical phenomena occurring in nonhnear systems which are 
too complex to analyze them directly, billiards have become popular model systems in many 
areas of physics: the main concepts of quantum chaos could be developed in the context of 
billiards {ll] (and refs. therein), they stimulated the foundations of modern semiclassics [l| 
and could be related to the Riemann hypothesis 121]. Even a justification of a probabilistic 
approach to statistical mechanics can be based on billiards Recently, they could 

be related to the physics of graphene J5| and have been exploited in order to improve the 
efficiency of thermoelectric materials |6|. Meanwhile, they are experimentally accessible in 
a variety of setups such as microwave resonators , |8[, atom-optical- and semi- 

conductor heterostructure based devices [3] as well as quantum dots 13|. 
In recent years particularly time-dependent billiards attracted much attention. A related 



key topic is the 
acceleration 14 



ong-standing quest of the general requirements of the occurrence of Fermi- 



16| (and refs. therein), i.e. under which conditions does a system 'heat up' 



unboundedly due to a driving force? The difficulty of a general answer becomes already 



clear by the exemplary findings, that correlations can either cause the 



acceleration, as in the one- dimensional Fermi-Ulam model 
exponential acceleration 19l-|22|. 
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suppression of Fermi- 



18| or can even yield 



Due to their conceptual simplicity and the complexity of their phenomenology, time- 
dependent billiards are perfectly suited to employ ideas which are based on characteristic 
phenomena of complex systems: In [23!] the fact that the set of attractors in the phase space 
of a dissipative nonlinear systems determines its possible asymptotic behavior was employed 
to an elliptical billiard with a time-dependent boundary in order to achieve a filter for 
neutral particles, working simultaneously in momentum and configuration space. Cutting a 
small hole into the boundary a synchronized particle beam is obtained. Such a dynamical 
particle filter, based on general concepts of nonlinear dynamics provides a simple, versatile 
and efficient alternative to the extensively used standard filtering-devices such as lenses, 
nozzles, mirrors or collimators. However, the corresponding results were achieved using 
a dissipative elliptical billiard with oscillating hard boundaries which are experimentally 
challenging to realize and turn out to lead to relatively high loss rates in terms of particles 
stopping inside the billiard. 

Static boundaries are comparatively simple to realize and design and lead in combination 
with a time-dependent scatterer in the interior of the billiard to an alternative dynamical 
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system. This raises the question, whether particle focusing as demonstrated in 23j can 



23| concerns the relatively high 



be obtained also in such a setup. A further limitation of 
particle losses due to rotators in the phase space of elliptical billiards [23 1. 
The main goal of this work is to provide and analyze a model which combines experimental 
simpleness with high-efficiency asymptotic focusing and synchronization of random ensem- 
bles of neutral particles on a single periodic trajectory. A billiard with oval boundaries 
and oscillating elliptical scatterer in the interior is shown to provide the appropriate setup 
and demonstrated to yield the desired focusing effects for sufficiently slow driving and an 
elongated shape of the inset. 

Motivation of the Model It is first discussed why for the class of billiards with a static 
boundary and a time-dependent scatterer in the interior, an oval for the former and an 
ellipse for the latter appears to be the simplest promising setup in order to achieve a highly 
efficient particle focusing. Then, a corresponding model is proposed. The obvious starting 
point would be to replace the breathing elliptical boundaries in [23] by static ones and to 
place a circular, oscillating scatterer in the interior of the billiard. However, due to the 
presence of rotators in the elliptical billiard j3] relatively high particle losses would be 
;he consequence which can be only avoided by introducing local defects at the boundary 



23| . Thus, the natural extension is to consider oval boundaries which in the following 
turn out to be perfectly tunable to avoid particle losses. Concerning the geometry of the 
inset, the natural starting point of a circular inset turns out to result either in attractors of 
chaotic type, i.e. in a non-synchronized, non-periodic asymptotic dynamics or in very low 
focusing efficiencies. These considerations therefore result in the following model: Consider 
the dynamics of non-interacting point particles of mass m in a two-dimensional static oval 
billiard with an oscillating hard scatterer of elliptical shape in the interior (Fig. [1]) in the 
presence of a velocity-proportional friction force = k ■ v. The static oval boundary is 



defined as a plain, closed curve xb(0) parametrized by < < 27r [25| : 



= [((5/2 + 1) sin0 + (5/6) sin(30)] 
ysicj)) = [{6/2 - 1) COS0 - (5/6) cos(30)] (1) 

Therein < 5 < 1 is a measure for the deviation of the oval from the circle (which is 
reobtained for 6 = 0) and for the maximal curvature, occurring at = 0, vr. The respec- 
tive boundary of the oscillating inset can be defined directly by the following coordinate 



FIG. 1. Sketch of the oval bilhard with the oscillating inset in its interior at two different collision- 
times (full and short dashed purple line), a typical trajectory (red, long dashed) and the angles 
a„ , /3n for an exemplary collision with the static boundary. 



representation 



[xi-m? , [yi-mr 




f^(t) = acos{uJ); fy{t) = bcos{uyt), 



(2) 



where A, B are the semi-major and semi-minor axes of the ellipse and u^, Uy, a, b denote the 
frequencies and amplitudes of the oscillation in x, y-direction, respectively. 
The corresponding free-space dynamics of a particle with coordinates x = {x,y), initial 
position xq and velocity vq is described by the Newtonian equation of motion mx + fcx = 0, 
with the solution 



Thus, particles inside the billiard move on straight lines between successive collisions with 
boundary /inset and their dynamics can be captured by a mapping from one collision to the 
next one. Collisions can of course happen in any order with the boundary and inset and not 
only in an alternating way. Particles which do not hit the driven inset for a sufficiently long 
time eventually stop and are considered as lost. 

In order to minimize losses, the occurrence of rotors in the phase space, corresponding to 
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x(t) = xo + ^ (1 - e-C^/-)*) 



(3) 



particles traveling around the inset without hitting it, must be suppressed. Accordingly, in 
the following a large deviation of the oval from the circle {6 = 0.99) and mainly elongated 
insets are chosen (cf. [T]). 

Mapping: A mapping which describes the dynamics of particles in the billiard from 
collision to collision is formulated in the following. Let n count the total number of collisions 
of the particle with the boundary /inset and )f:„,x„, v„ denote time, position and velocity at 
the instance after the corresponding collision. Then, the dynamics of the system can be 
described by a five-dimensional, non-area-preserving and non-Hamiltonian mapping 

M : (t„,x„,v„) ^ (tn+i,Xn+i, Vn+i) (4) 

Since it is not clear a priori whether the respective next collision happens with the boundary 
or with the inset, both are calculated and the smaller one is selected. In the following the 
iteration of the mapping is discussed for the possible cases, collisions with the oval boundary 
and the elliptical inset. 

Collision with the oval boundary: First on, the slope of the trajectory after the respective 
last collision is calculated, then position and velocity and finally the corresponding time of 
the next collision follow. Let a„ measure the angle between the tangent to the boundary at 
the n-th collision point and the and /3„ the angle between the velocity vector and 

the X-axis at this collision point (cf. [1]). Then, the slope of the tangent at the n-th collision 
point mB,n can be calculated by 

dys/d^ 



mB,n 



da;B(0) 



dxB/d(j) 



;i -5/2)sin0+ (5/2)sin(30) 

(5) 



(1 + 5/2) cos + (6/2) cos(30) 

and an = arctan(mB_„) holds. With = 2a„ — /3„ the trajectory after the n-th collision 
has the slope m^.n := tan(/3„+i), i.e. it evolves along a straight line with a known slope of 
y{x) = mT,n{x — Xn) + Vn- Inserting the boundary-functions Xb{4>) and yB{(p) (Eq. [I]) for x 
and y yields an implicit equation determining the value of the angle at which the particle 
hits the boundary (in the inset-free) case at the respective next collision: 

g{(t)) : = mT,n [{S/2 + 1) sin0 + (5/6) sin(30)] 

- [{6/2 - 1) cos - {6/6) cos(30)] + (y„ - mT,nXn) = (6) 

This implicit equation determines the boundary parameter 0„+i = corresponding to the 
respective next collision point of particle and oval and can be solved numerically with a 



standard Newtonian root finding method. Inserting 0n+i into Eq. ([T]) yields the coUision 
points in configuration space x^_^i = XB{4>n+i) and Vn+i = XB{(pn+i)- Therein the superscript 
B denotes the colhsion with the static boundary. With Eq. Oand x„_|_i — x„ being parallel to 
the corresponding collision time results as t^^^ = {—m/k) In [1 — A;|x„+i — x„|/(mf„)] + 
tn- The new velocity v^+i directly follows with v^^-^ = f„ ■ exp [(— A;/m)(t„+i — tn)] and 
t^^+i/t^^+i = tan(/3"+i) = mT,n- 

Collision with the oscillating elliptical inset: In contrast to the oval boundary, for calculating 
the collision time with the inset, advantage can be taken from the fact, that a coordinate 
representation for the ellipse exists. Inserting Eq. with vq i— ?■ {v^^n, Vy^n) and 1 1— tn+i —tn 
component-wise into Eq. [2] yields the following implicit equation determining the time tn+i 
of the respective following collision with the inset: 



xn + ^(1 - e-^(W-*")) - f.(t,+i) ■ /A' 
k J 

= : = 0; - t„ > (7) 

Note, that Eq. (JTj) only has a finite solution when /(t„,+i) changes sign at least once in the 
interval (t„,t^+i). In this case, the respective smallest solution determines tn+i and the 
respective next colhsion occurs with the inset. Otherwise tn+i = x,„_|_i = x^.^^^, v„+i = 
v^_^j^ is chosen. For a collision with the inset, the velocity changes as 

v„+i = v„e[(-'^/™)(*"+^-*")l - 2 [n ■ (v„e[(-^/™)(*"+^-*")l - 5vi{t))] ■ n (8) 

where S\j(t) = ^tyiI(t)\^^^^^^_^ is the velocity vector of the inset at the respective collision 
point and n = (5v/(t)/|5v/(t)| is the corresponding normal vector. After each collision the 
slope can be directly updated via mT,n+i = "y^+i/'W^+i such that the respective part can be 
skipped above after a collision with the inset. All considerations concerning the time of the 
respective next collision of course hold independently of the sequence (boundary/inset) in 
which the collisions occur and the mapping can be straightforwardly iterated. 

Particle Focusing The mapping is numerically iterated for ensembles of 10^ initial con- 
ditions for 2.8 ■ 10^ collisions and various choices of inset geometry and driving law. For 
the following discussion four representative examples are chosen. Initial velocities are cho- 
sen randomly with uniform probability in f G [5, 20] (the lower offset is to avoid particles 
stopping before hitting the inset for the first time) . Regarding the evolution of the ensemble 
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FIG. 2. Time-evolution of the ensemble-averaged particle velocity. Parameters uj^ 



1.0: 



Circle: A = B = 0.5; a = b = 0.1; k = 0.002; Ellipse 1: A = 1.0; B = 0.25; a = 0; b = 0.25; k 
0.002; Ellipse 2: a = 1.3; 6 = 0.05; A = 0;B = 0.05; k = 0.001; Ellipse 3: a = 1.2; 6 = 0.05; yl 
0.1; B = 0.05; k = 0.001. 



averaged magnitude of the velocity < >, according to Fig. [2], saturation is observed for 
all four visualized combinations of geometry and driving law. This saturation suggests that 
all particles have approached attractors, i.e. that the system is in its asymptotic state at 
the end of the simulation-time. For the circular inset (circle) and the relatively fast driven 
elliptical inset (ellipse 1) the average-velocity grows initially, then saturates at a certain 
value well above the average initial velocity. In contrast to that the average velocity de- 
creases for both slowly driven elliptical billiards (ellipse 2,3) and converges towards values 
of < \vn\ >~ 0.4. Comparatively strong driving (for example A = B = 0.3 for the circular 
inset - not shown) leads to unbounded growth of energy on the examined timescales and is 
thus inappropriate in order to achieve particle focusing. 

Now the asymptotic dynamics is analyzed in more detail. Fig. [3] (Circle) shows that the 
slowly driven circular inset leads to a broad asymptotic velocity distribution. Thus, it is 
inappropriate in order to achieve a focusing of particles on single trajectories. The broad dis- 
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FIG. 3. Velocity-distributions after 2.8 • 10^ collisions. Parameters like in Fig. [2j 
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FIG. 4. Time-dependent magnitude of the velocity after 2.8 • 10^ collisions for the configuration of 
ellipse 2 for two different initial conditions. Parameters like in Fig. [2j 



tribution together with the observed saturation of the mean ensemble velocity is indicative 
of the presence of a strange attractor, i.e. that trajectories are asymptotically focused onto 
a submanifold with lower (and non integer) dimension than the phase space. An obvious 
conclusion would be to assume that circular insets do not allow for the presence of isolated 
attractive periodic orbits since defocussing is stronger than the focusing of the oval bound- 
ary in regions where attractors are conceivable. Thus, elliptical insets of elongated shape 
should be more promising. However, the velocity distribution (Fig. |3l ellipse 1) shows that 
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even the combination of a very flat elliptical inset and moderate driving is not sufficient in 
order to achieve particle focusing: again a stationary, i.e. non-spreading, but broad velocity 
distribution and corresponding chaotic motion of the respective trajectories is obtained. 
A slowly and purely vertically driven flat elliptical inset provides an asymptotic velocity 
profile with two clearly pronounced peaks at \v\ ~ 0.2 and \v\ ~ 0.4 containing 92.7% of all 
initial conditions (Fig. [3l ellipse 2). Less than one percent of all particles are lost (low peak 
at |f I = 0). The peak at |f | ~ 0.2 corresponds to a periodic orbit attractor of period two. 
The corresponding asymptotic time-evolution of the magnitude of the velocity is shown in 
Fig. m (upper plot). The amplitude of the oscillations of \v{t)\ continuously decreases as 
the attractor is approached. The related motion in configuration space is a two-periodic 
bouncing exactly on the y-axis between inset and boundary. About 39.3% of all initial con- 
ditions are asymptotically focused on this periodic orbit and are thereby synchronized to 
perform asymptotically one and the same dynamics. The peak at |f | ~ 0.4 corresponds to 
a periodic orbit attractor of period four (Fig. HI lower plot) which splits into four individual 
peaks in a very fine resolved version of Fig. [3] (ellipse 2) which merge in the long-time limit. 
About 53.4% of all initial conditions are focused on this periodic trajectory and are thereby 
perfectly synchronized in the long-time limit. In configuration space, the corresponding 
particles are focused on the same orbit as for the |f | ~ 0.2 peak. The remaining particles 
are focused on a variety of more complicated periodic trajectories, all in the velocity regime 
0.2 < \v\ < 0.39 and very few ones on a trajectory around \v\ ~ 5 (Inset in FiglH ellipse 2.). 
This demonstrates, that it is indeed possible to achieve the focusing of randomly distributed 
ensembles of neutral particles on periodic trajectories with close to perfect efficiency by us- 
ing static oval billiards with elongated, slowly driven elliptical insets. 

It is now suggestive that the respective, purely vertically aligned periodic orbit can also 
become attractive for a circular inset but with a much smaller basin of attraction. Since the 
circle is stronger defocussing than the boundary is focusing sophisticated time-correlations 
are required in order to approach it. Corresponding numerical simulations confirmed this 
expectation: for sufficiently slow driving, circular insets also yield attractive periodic tra- 
jectories and may thus be used in principle for particle focusing, but the corresponding 
efficiencies are extremely poor. The best achieved efficiency was well below 10%, i.e. typi- 
cally most particles stop before reaching periodic orbits. 

Further numerical simulations have been performed for a variety of combinations of slow 
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driving laws and elongated elliptical insets and always yielded both, the occurrence of fixed 
point attractors in phase space and small loss rates. Consequently the demonstrated particle 
focusing does not require any fine tuning of parameters. Fig. [3] (Ellipse 3) exemplarily shows 
the resulting velocity distribution for a slowly driven, vertically oscillating inset. Again, 
more than half of the ensemble is focused asymptotically onto the above-discussed period- 
four attractor, but the rest of the ensemble asymptotically follows a variety of periodic orbits 
with velocities between |f | ~ 0.18 and |f | ~ 0.55. A very different velocity profile could be 
achieved for example for a = 1.0; b = 0.05; A = 0.3; B = 0.0125; k = 0.001; = 2.5; Uy = 0.5 
but only for the price that a relevant part of the ensemble is asymptotically governed by a 
chaotic attractor and thus remains non-synchronized. 

Note that the particle-focusing mechanism is governed solely by the topological structure of 
the phase space and thus works independently of the specific choice of the initial conditions. 
Key ingredients for high efficiency focusing are strong curvature of the oval boundary at 
certain positions and elongated insets such that particles are trapped in the upper or in the 
lower half of the billiard for many collisions. 

Finally, possible experimental implementations are discussed. The presented numerical re- 
sults should be accessible to particle-clusters in a micromechanical realization of the billiard- 
system with medium. Alternatively, a purely optical verification of the mechanism should 
be equally possible by using atoms in laser-based optical traps of oval shape, counter- 
propagating laser fields yielding friction due to the optical molasses and an oscillating blue 
detuned laser taking the role of the oscillating, reflecting inset. According to the following 
consideration the parameter sets used in this work should be effectively accessible in both 
setups: In Eqs.(|3]),([6]),([7]),([8]) the length scale occurs only in the driving amplitudes a, b and 
in aux,buy, the time-unit appears in the driving frequencies Ux,ujy, in aux,buy and in the 
friction coefficient, while the mass only appears reciprocal to the friction coefficient which 
scales with mass and time scale. A desired friction coefficient can be effectively realized by 
choosing particles with appropriate mass or also via the frequency of the driving, but the 
latter determines the size of the billiard, which otherwise can be chosen arbitrarily, as long 
as aUx, buy are preserved. 

Conclusion It has been demonstrated, that dynamical focusing of ensembles of neutral 
particles in energy and configuration space, is possible in billiards with static boundaries 
and an oscillating scatterer in the interior. Employing appropriate oval boundaries and 
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sufficiently slowly driven elongated elliptical insets, close to perfect efficiency can be achieved. 
The suggested focusing device can be implemented either in a micromechanical or in a 
purely optical setup yielding a versatile alternative to common focusing devices like lenses 
or collimators. 
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